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Abstract 

We extend the results of a previous analysis of ours showing that, when both heavy 
and light flavour effects are taken into account, successful minimal (type I + ther- 
mal) leptogenesis with 5'O(10)-inspired relations is possible. Barring fine tuned 
choices of the parameters, these relations enforce a hierarchical RH neutrino mass 
spectrum that results into a final asymmetry dominantly produced by the next- 
to-lightest RH neutrino decays dominated leptogenesis). We present the con- 
straints on the whole set of low energy neutrino parameters. Allowing a small 
misalignment between the Dirac basis and the charged lepton basis as in the quark 
sector, the allowed regions enlarge and the lower bound on the reheating tempera- 
ture gets relaxed to values as low as ~ 10^'' GeV. It is confirmed that for normal 
ordering (NO) there are two allowed ranges of values for the lightest neutrino mass: 
mi ~ (1 - 5) X 10"^ eV and mi ~ (0.03 - 0.1) eV. For mi < 0.01 eV the allowed 
region in the plane ^13-^23 is approximately given by 6*23 < 49° + 0.65(6*13 — 5°), 
while the neutrinoless double beta decay effective neutrino mass falls in the range 
mee = (1 - 3) X 10-3 eV for 6*13 = (6° - 11.5°). For mi > 0.01 eV, one has quite 
sharply m^e ~ rni and an upper bound ^23 ^ 46°. These constraints will be tested 
by low energy neutrino experiments during next years. We also find that inverted 
ordering (10), though quite strongly constrained, is not completely ruled out. In 
particular, we find approximately ^23 — 43° -|- 12° log(0.2 eV/mi), that will be fully 
tested by future experiments. 



1 Introduction 



With the discovery of neutrino masses and mixing in neutrino oscillation experiments, 
leptogenesis [H [2] has become the most attractive model of baryogenesis to explain the 
observed matter-antimatter asymmetry of the Universe. This can be expressed for exam- 
ple in terms of the baryon-to-photon number ratio and is very well measured by CMB 
observations [3] to be 

r/g^^ = (6.2±0.15) X 10-^°. (1) 

Leptogenesis originates from the see-saw mechanism [3] that is based on a simple extension 
of the Standard Model where right-handed (RH) neutrinos with a Majorana mass matrix 
M and Yukawa couplings h to leptons and Higgs are added. Within S'0( 10) models, three 
RH neutrinos Ni {i = 1,2,3) are nicely predicted and for this reason they are traditionally 
regarded as the most appealing theoretical framework to embed the seesaw mechanism. 

However, within the simplest set of assumptions inspired by 5'O(10) models [5], barring 
strong fine-tuned degeneracies in the RH neutrino mass spectrum and using the experi- 
mental information from neutrino oscillation experiments, the traditional A'^i-dominated 
leptogenesis scenario predicts an asymmetry that falls many orders of magnitudes below 
the observed one ^ i7j. This is because, within A^^i -dominated leptogenesis, where the 
spectrum of RH neutrinos is hierarchical and the asymmetry is produced from the decays 
of the lightest ones, successful leptogenesis implies a stringent lower bound on their mass 
[8], Ml > O(IO^) GeV. On the other hand, 5*0(10) grand-unified theories typically yield, 
in their simplest version and for the measured values of the neutrino mixing parameters, 
a hierarchical spectrum with the RH neutrino masses proportional to the squares of the 
up-quark masses, leading to Mi = O{10^) GeV and therefore to a final asymmetry much 
below the observed one. 

However, it has been shown [9] that, when the production from the next-to-lightest RH 
neutrinos [10] and lepton flavour effects [11] are simultaneously taken into account |12j . 
the final asymmetry can be generated by the decays of the next-to-lightest RH neutrinos 
and allowed regions in the low energy neutrino parameter space open up. 

In this paper we proceed with the analysis of [9j and present the resulting constraints 
on all low energy neutrino parameters. The paper is organized as follows. In Section 2 
we discuss the current experimental status on low energy neutrino parameters, we set up 
the notation and describe the general procedure to calculate the the asymmetry and find 
the constraints. In Section 3 we first consider the case already studied in [9], when the 
Dirac basis and the charged lepton basis coincide and then, in Section 4, we allow for 
a misalignment between the two bases not larger than that one described by the CKM 
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matrix in the quark sector. Finally, in Section 5 we present a global scan in the space of 
parameters where all possible cases between the case of no misalignment and the case of a 
misalignment at the level of the CKM matrix are taken into account. We also discuss two 
scenarios, one at small mi and one at large mi, and show how, within 5'O(10)-inspired 
models, minimal leptogenesis could be tested in future low energy neutrino experiments. 

Notice that our discussion is made within a non-supersymmetric framework. Recently 
a study of S'O(10)-inspired models within a supersymmetric framework has also enlight- 
ened interesting potential connections with lepton flavour violating decays and Dark Mat- 
ter [13]. An analysis of leptogenesis within left-right symmetric models, where a type II 
seesaw contribution to the neutrino mass matrix is also present, has been performed 
in [H]. Within these models, the minimal type I scenario considered here represents a 
particular case recovered under specific conditions. 



2 Experimental information and general setup 

After spontaneous symmetry breaking, a Dirac mass term itld = hv, is generated by the 
vacuum expectation value (VEV) t> = 174 GeV of the Higgs boson. In the see-saw limit, 
M ^ m^i, the spectrum of neutrino mass eigenstates splits in two sets: three very heavy 
neutrinos, Ni,N2 and respectively with masses Mi < M2 < M3 almost coinciding 
with the eigenvalues of M, and three light neutrinos with masses mi < m2 < m^, the 
eigenvalues of the light neutrino mass matrix given by the see-saw formula ^ , 

It 

m^ = -m£,— — m^, (2) 

that we wrote in a basis where the Majorana mass matrix is diagonal defining Dm = 
diag(Mi, M2, M3). The symmetric light neutrino mass matrix m^ is diagonalized by a 
unitary matrix [/, 

f/t W = (3) 

with Dm = diag(mi, m2, ma), that, in the basis where the charged lepton mass matrix is 
diagonal, can be identified with the lepton mixing matrix. 

Neutrino oscillation experiments measure two neutrino mass-squared differences. For 
NO one has m| — m| = Am^^^ and m| — m^ = Am^^^j. The two heavier neutrino masses 
can therefore be expressed in terms of the lightest neutrino mass mi as 

m2 = ^mj + m^^i , and mg = ymfTm^^, (4) 

where we defined matm = \/ ^'^atm + ^^Li = (0.050 ± 0.001) eV and m^oi = \/ Am^^j = 
(0.00875 ±0.00012) eV Recently, a conservative upper bound on the sum of neutrino 
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masses, < 0.58 eV (95% CL), has been obtained by the WMAP collaboration [3] 

combining WMAP 7 years data plus baryon acoustic oscillations observations and the 
latest HST measurement of Hq. Considering that it falls in the quasi-degenerate regime, 
it straightforwardly translates into 



mi < 0.19 eV (95% CL) . 



(5) 



We will adopt the following parametrization for the matrix U in terms of the mixing 
angles, the Dirac phase 6 and the Majorana phases p and cr [16] 



U 



Cl2 Ci3 



"Sl2 Ci3 



-Sl2 C23 - Ci2 S23 S13 c'^ C12 C23 - S12 S23 S13 c"^ 



Sue-'' \ 

S23 Cl3 



diag ( 



\ S12 S23 - C12 C23 Si3 e"^ -C12 S23 - S12 C23 Si3 e' ' C23 Ci3 / 



(6) 



and the following 2 a ranges for the three mixing angles [15 



= (31.3° -36.3°), 023 = (38.5° - 52.5°) , ^13 = (0° - 11.5°) . (7) 
In the case of 10 the expression of m2 in terms of mi becomes 



m2 = x/mi + ml^ 



m. 



sol ' 



(8) 



while the expression for m3 does not change. With the adopted convention for the light 
neutrino masses, mi < m2 < m3, the case of 10 corresponds to relabel the column of the 
leptonic mixing matrix performing a column cyclic permutation, explicitly 



U 



-,-i5 



C12 Ci3 



/ si3e 

S23 Cl3 

\ C23 Ci3 S12 S23 — C12 C23 Si3 e 



S12 C23 ~ C12 S23 Si3 e 

i<5 



i(5 



S12 Ci3 
C12 C23 ~" Si2 S23 Si3 e 



\ 



i(5 



"C12 S23 — S12 C23 Si3 & ' J 



diag (e^'",e''',l) . 

(9) 

The predicted baryon-to-photon ratio rjB is related to the value of the final {B — L) 
asymmetry N^b-l by [2] 

riB - 0.96 X 10-2iV£_^ , (10) 

where Nb-l is the B — L number in a co-moving volume that contains on average one 
RH neutrino Ni in thermal ultra-relativistic equilibrium abundance (T ^ Mi). 
The Dirac mass matrix can be diagonalized by a bi-unitary transformation 



rriD = Vl D.mo Un , 



:ii^ 
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where -Dm^, = diag(Ai, A2, A3). The matrix Ur can be obtained from Vl, U and m^, 
considering that it provides a Takagi factorization [17] of fl8\ |5] 

^ D7^\ Vl U Vl D~l , (12) 

exphcitly 

M-^ = UnD-^Ul. (13) 

For non degenerate Mj, the matrix Vr can be determined noticing that it diagonahzes 
M-i (M-^)t, i.e. 

M-\M-y = UnDl^Ui. (14) 

This relation determines [/r unless a diagonal unitary transformation, since any Vr = 
UrD~^^ is also a solution. However, given a Ur, one can fix from the eq. f|T3|) . 

= ^JDMU^M-^U*^ (15) 

and in doing so ^7/^ is unambiguously determined. Inspired by SO{10) relations, we can 
parameterize the eigenvalues of rriD in terms of the up quark masses as 

Ai = aim„, X2 = a2mc, X3 = asmt. (16) 

Within SO{10) models one can expect ai = 0{1) and we will refer to this case. The reader 
is invited to read Ref. [9] for a more comprehensive discussion about these S'O(10)-inspired 
relations. Notice however that our results will be valid for a much broader range of values, 
since, quite importantly, it turns out that they are independent of ai and 0:3 provided 
M3 ^ M2 and Ml < 10^ GeV. With the parametrization eq. (JT6II and barring very special 
choices of parameters where the RH neutrino masses can become degenerate ^ □, the RH 
neutrino mass spectrum is hierarchical and of the form (for generic expressions in terms 
of the low energy parameters, see Ref. [7j) 

Ml : M2 : M3 = (ai mj^ ; (^2 m,)^ : (^3 • (17) 

As we said, the values of ai and are actually irrelevant for the determination of the 
final asymmetry (unless ai is unrealistically large to push Mi from ~ 10^ GeV above the 
lower bound ~ 10^ GeV to achieve successful A''i leptogenesis) . On the other hand, the 

^ As in 1^, we consider only solutions where M3/M2 and M2/M1 > 10. This is clearly a conservative 
condition, since the asymmetry gets enhanced when M2 — M3 or M2 — Mi. However, in this way, 
we only neglect very special points in the parameter space yielding M3/M2 and M2/M1 < 10. We will 
comment again later on this point. 
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value of a2 is relevant to set the scale of the mass M2 — 2(a2^c)^/'^3 (valid for ^13 ^ 0) 
of the next-to-lightest RH neutrino mass, but it does not alter other quantities crucial for 
thermal leptogenesis, such as the amount of wash-out from the lightest RH neutrinos. 
Defining the flavoured CP asymmetries as 

£2a = =^ , 18 

r2 + r2 ^ ' 

these can be calculated using [19] 



[1 + x) In 



(20) 



16n{hmh2 [ ' '''' 3(X3/X2 

(19) 

where 

■l-ha;\ 2-x 
X J 1 — X 

and is the decay rate of the RH neutrino N2 into the flavor a with couplings given 
by the Yukawa's matrix h. We will assume an initial vanishing A^2-abundance instead of 
an initial thermal abundance as in [9]. In this way, a comparison of the results in the two 
analyses gives a useful information about the dependence of the final asymmetry on the 
initial N2 abundance when successful leptogenesis is imposed. 
Let us now define the flavored decay parameters as 



where H is the Hubble rate 



|2 



= ~ 1.08 X 10"^ eV , 22 

3x/5 Mpi ^ ^ 

(7* is the number of the effective relativistic degrees of freedom and Mpi is the Planck 
mass. The total decay parameters are then just simply given by Ki = Kia- It is also 
convenient to introduce the quantities = -f^2a/-^2- 

From the decay parameters one can then calculate the efficiency factors that are the 
second needed ingredient, together with the CP asymmetries, for the calculation of the 
final asymmetry. These can be well approximated by the following analytical expression 
[20]@ 

n{K2, K2a) = nl{K2, K2a) + «:Ui^2, K2a) , (23) 



^It is in quite a good agreement with the numerical resuhs shown in |21) . The maximum difFerence is 
30% at the peak for ^ 1- For S> 1, the difference is below 10%. 
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where the negative and the positive contributions are respectively approximately given 
by 

kI{K,,K,^) ^ -A (^e1-A^^.(-.) _ 1^ , (24) 

N^,, (4^) ^ N{K) ^ . ^^^'^ , (25) 



and 



where 



1 + v/iV(K^ 



2.5 



^B(i^2a) ^ 2 + 4ir°a'e = 0(1 ^ 10) . (27) 

The 5'O(10)-inspired conditions = 0{1), yield a RH neutrino mass spectrum with 
Ml <^ 10^ GeV < M2 < 10^^ GeV ^ M3, though, as we already noticed, this spectrum is 
obtained for a broader range of values. In this situation, the asymmetry is dominantly 
produced from N2 decays at T ~ M2 in a two flavour regime, i.e. when final lepton 
states can be described as an incoherent mixture of a tauon component and of coherent 
superposition of a an electron and a muon component. Therefore, at the freeze-out of 
the N2 wash-out processes, the produced asymmetry can be calculated as the sum of two 
contributions, 

iV^rf 2 ~ e2r n{K2, K2r) + £2e+M '«(^2, K2e+^) , (28) 

where e2e+^ stands for ege+M = £2e + ^2^ and K2e+f, = K2e + K2^- 

More precisely, notice that each flavour contribution to the asymmetry is produced 
in an interval of temperatures between M2/[zb(-^2q) — 2] and M2/{zB{K2a) + 2], with 
a = r, e + /i. 

At T < 10^ GeV the coherence of the e + ^ quantum states breaks down and a three 
flavour regime holds, with the lepton quantum states given by an incoherent mixture of 
e, /i and r flavours. The asymmetry has then to be calculated at the A^i wash-out stage 
as a sum of three flavoured contributions. 

The assumption of an initial vanishing A^2-abundance allows to neglect the phantom 
terms in the muon and in the electron components [22] so that the final asymmetry can 
be calculated using the expression 

N's^L ^Le2e+M't(i^2e+M)e-^'^^'= + 7§^£2e+M't(i^2e+M) e"^ ^^-+£2r /t(i^2r) e"^ . 

(29) 
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Notice that successful leptogenesis relies on points in the parameter space where one 
out of the three Ki^ < 1. From this point of view the constraints on low energy neutrino 
experiments that we will obtain should be quite stable against effects that could enhance 
the asymmetry such as a resonant enhancement for special points where {M^ — M2) / M2 ^ 
1. Such effects are however still able to relax the lower bound on M2 and on the Trh? 
since the ilia's do not depend on M2. 

3 The case Vl = I 

We start from the case Vl = I that has been studied already in [9] deriving constraints in 
the plane mi —^13 for NO. Here we show constraints on all low energy neutrino parameters, 
including the case of 10. 

3.1 Normal ordering 

Let us first discuss the case of NO. In Fig. 1 we plotted the final asymmetry rjB for the 
same three sets of values of the involved parameters as in the Fig. 4 of Ref. [9], where 
these three choices were corresponding to three different kinds of solutions for successful 
leptogenesis. This time the third solution (right panel), is suppressed and successful 
leptogenesis is not attained. In [9], this was the only solution corresponding to a final 
asymmetry dominantly in the muon flavour instead than in the tauon flavour (as for the 
first two). The suppression that we find now is explained partly because we are adopting 
a correct determination of the phases in the Ur matrix (cf. eq. f|T5|) ) and partly because 
we are now assuming an initial vanishing Ai"2-abundance instead than an initial thermal 
one. We will see however that, allowing for Vl 7^ /, this kind of solution will again yield 
successful leptogenesis in some allowed regions of the parameter space, characterized in 
particular by large values of mi ~ 0.1 eV. 

The solution in the central panel is also partly suppressed and successful leptogenesis 
is not attained. However, in a parameter scan, we find that this kind of solution can 
still give successful leptogenesis for slightly different values of the parameters than those 
indicated in the figure caption. In this case the difference with respect to the results in 
[9] is explained just in terms of the different assumption on the initial abundance. This 
dependence on the initial conditions is due to the fact that ~ 1, i.e. the solution falls 
in the weak wash-out regime at the production. 

Finally, the first solution (left panel) is fully unchanged. It therefore exhibits a full 
independence of the initial conditions and this is in agreement with the fact that it respects 
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Figure 1: Case Vl = /, NO. Plots of the relevant quantities for three choices of the involved 
parameters as in Fig. 4 of Ref. [9J: ^^13 = 5°, ^23 = 40° 612 = 33.5° in all three cases. The 
values of the phases are different in the three panels (radiants): 6 = a = 0, p = 1.5 (left); 
S = 5.86, p = 0" = 3 (center); 6 = vr/3,p = 0.02,0" = 7i/2 (right). The long-dashed red 
lines correspond to a = r, the dashed blu§ lines to a = /i and the short-dashed dark 
yellow lines to a = e. In the bottom panels the horizontal dotted line is the 2a lowest 
value T]'^^^ = 5.9 x 10~^° (cf. ([1])), the solid line is 77^ while the dashed line is r/^'"*^^. 



all the necessary conditions for the independence on the initial conditions found in 
Notice that these conditions also enforce an efficient wash-out of a possible pre-existing 
asymmetry. 

A scan in the space of parameters confirms that these three solutions obtained for 
special sets of values are actually representative of the three general kinds of solutions 
that come out and, therefore, the drawn conclusions apply in general. 

In the panels of Figure 2 we show the results of such a scan that highlight the allowed 
regions in the parameter space projected on different two-parameter planes. The scatter 
plots have been obtained scanning the three mixing angles ^12,^23 and ^13 over the 2a 
ranges eqs.([7]), the three phases S^p^a over the ranges [0, 27r] and the absolute neutrino 
mass scale for mi < 1 eV. These ranges also coincide with those shown in the plots, 
except for mi where the plots are for mi > 10~^ eV simply because no allowed solutions 
have been found for lower values. The shown results have been obtained in two steps. A 
first scan of 0(10^) points has yielded a first determination of the allowed regions. With 
a second scan of additional 0(5 x 10^) points, restricted to the excluded regions, we have 
then more robustly and sharply determined the contours of the allowed regions. Notice 
that these regions have no statistical significance and the random values of the parameters 
have been generated uniformly. 

In the top left panel the three RH neutrino masses are plotted versus mi. We have 
also plotted the lower bound on the reheating temperature calculated as 

T™" ~ — \ . (30) 

This calculation relies on the fact that in the case Vl = I the solutions, as we commented, 
always fall in a tauon A^2-dominated scenario. It can be seen that the lowest bound is 
given by ~ 2 x 10^° GeV that in a supersymmetric version, if unchanged, would be 
marginally reconcilable with the upper bound from the gravitino problem [25- This is 
another reason to extend our investigation to cases with V/^ 7^ / in next sections. 

In the top central panel we have then plotted the allowed region in the mi — 9i^ plane 
that can be compared with an analogous figure in [9]. Here, however, we show only those 
points that respect the condition rjB > 5.9 x 10^^° but for 2 different values of 0:2 = 4, 5 
(in [9] we were only showing points for 0:2 = 5). The (red) star represents a point found 
for a minimum value 02 = 3.4. This point basically roughly indicates where the maximum 
of the asymmetry occurs in the parameter space for a fixed value of a2- We will continue 
to use this convention (yellow circles for 0:2 = 5, green squares for 02 = 4 and red stars 
for minimum found 0:2 value) throughout the next figures. The structure of the allowed 
region in the mi — 613 plane can be understood as follows. Since oc (M2/M3) and 
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Figure 2: Case Vl = /, NO. Scatter plot of points in the parameter space that satisfy the 
condition rj^ > 5.9 x 10^^ for three values of the crucial parameter a2'- 0^2 = 5 (yellow 
circles); a2 = 4 (green circles); 02 = 3.4 (red stars). In the top left panel the lower bound 
on TptH (cf. eq. ([30])) is also indicated for the same values of 0:2 but with different symbols: 
«2 = 5 (grey squares), ^2 = 4 (black squares), 02 = 3.4 (blue star). The three mixing 
angles are in degrees, the three phases in radiants. The dashed line in the left central 
panel is the eq. 
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Ms oc m^^, we immediately deduce that a large lepton asymmetry in the tau flavor may 
be produced only for sufliciently large values of mi. This is rather easy to understand. 
If mi tends to zero, we go into the so-called decoupling limit, M2/M3 ~ 0. As the CP 
asymmetry needs (at least) two heavy states to be generated at the one-loop level, and 
disregarding the contribution from the A^^i, must vanish. The wash-out parameter 
is (9(25) [9j and therefore the flnal baryon asymmetry may be estimated to be 

r]B^5x 10-3 £2. ^ 5 f ■ 10-10 , (31) 



mi> { — ] 10-3 eV, (32) 



which requires 

rrii > I 

for NO. This estimate holds if the wash-out from the interaction with A^^i is negligible, i.e. 
Kir ~ 1- Of course, the smaller is mi, the smaller Ki^ needs to be. For mi = O (10-^) 
eV, the only possibility is that Kir is signiflcantly below unity. Extending the analysis of 
Ref. one flnds 

o A ^ m2tan6'23 „ „ . .„„^ 
Si3 cos (5 - 2cr) > -= ~ 0.04. (33) 

3v 2m3 

To get the feeling of the flgures involved, we may set 5 ~ 2cr and flnd that the wash-out 
mediated by the A'^i's vanishes for an experimentally allowed value of the mixing between 
the flrst and the third generation of LH neutrinos, ^^13 > 2.3° in agreement with our 
numerical results. If mi is larger than O (10-^) eV, then Kir = C^(l) is allowed and ^13 
can be taken to be vanishing. Notice also that the lower bound eq. (133|) on 6*13 increases 
with tan^23- This nicely reproduces the linear dependence emerging from the numerical 
results in the left column middle panel for the plane ^13 — ^23 and that is described, roughly 
for 0^2 = 5 and more accurately for 02 = 4, by 

^^23^ 44° + 4 (^13 -7°), (34) 

represented with a dashed line in the panel. In the top right panel we show the allowed 
region in the plane mi — 6*23. 

The CP non conserving terms in neutrino oscillation probabilities can be expressed in 
terms of the Jarlskog invariant Jcp given by [25] 

Jcp = Im[f/^3f/e2f/;2f/:3] (35) 
= C12 S12 C23 S23 Ci3 Si3 sin 5, (36) 

such that 

, ^ fAml-L\ 
Pu^^ua - Pu^^up =4:Jcp Sap;kj sin ( 1 , (37) 

k>3 ^ ^ 
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Figure 3: Case Vl = I, NO, 0^2 = 4. Constraints on the mixing angles obtained without 
imposing the current experimental information from neutrino oscillation experiments (blue 
points) compared to those previously obtained (green points). Notice that the regions 
exhibit a tt periodicity and they are specular around 7r/2 so that all mixing angles can be 
limited to the physical range [0, 7r/2]. This can be proven to hold on very general grounds 
[25] and therefore this plot can be regarded as a consistency check as well. 

where Sap-kj = ±1- In the bottom left panel we show the allowed points in the plane 
Jcp — diz- It can be noticed that a non zero value of Jcp is not crucial. Looking at the 
bottom-central panel, it is interesting to notice that the allowed regions for the Majorana 
phases are centered approximately around cr = nvr and p = {n + 1/2)tt. 

These play a role in the determination of the effective Majorana mass of in /3/30z/ 
decays that is given by 



{31 



= Imici^ci^e^''' + 7712 CI3 + 7713 e^'^^-^'l . (39) 

In the bottom-right panel one can see how there is a precise relation between and mi , 
given approximately by rngg — m-i. It can be also noticed that there is quite a strict lower 
bound TTiee > 1.5 X 10^^ eV. Lowest values ^ 2 x 10~^ eV are the most favoured ones 
in this case. Though current planned experiments will not be able to test the full allowed 
range, it is still interesting that they will test it partially, tightening the constraints on 
the other parameters as well. 

We have also made an interesting exercise. We determined the constraints without 
making use of any experimental information on the mixing angles and letting them just 
simply variate between 0° and 360°. The results are shown in Fig. [31 First, notice that the 
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lower bound on mi relaxes of a few orders of magnitude (see left panel). Then notice quite 
interestingly that small values ^13 < 10° are well allowed for mi > 10~^eV but values 
30° > ^13 ^10° would have been very marginally consistent. Therefore, the current 
bound 6*13 < 10° seems to match quite well with successful S'O(10)-inspired leptogenesis 
On the other hand, values 62^ ^ 30° would have been more optimal for ^13 < 10° than 
the current experimental large atmospheric values (see the central panel in the figure). 
However, they are still allowed thanks to the observed range of values of the solar neutrino 
mixing angle (see the right panel). For the solar neutrino mixing angle there is no real 
favourite range of values for ^13 < 10°. 



3.2 Inverted ordering 

Let us now discuss the results for 10. It has been shown |26] that in grand unified models 
with conventional type I seesaw mechanism one can always find, for any NO model satis- 
fying the low energy neutrino experimental constraints, a corresponding 10 model. There- 
fore, though they exhibit some unattractive features that quite strongly disfavour them 
(e.g. instability under radiative corrections), 10 models within grand unified theories are 
not unequivocally excluded. It is therefore legitimate to check whether the requirement of 
successful leptogenesis can somehow provide some completely independent information. 

We repeated the same scan performed in the case of NO and the results are shown 
in figure IH the analogous of the figure [2] for the NO case. One can see that 10 is only 
very marginally allowed. For 02 = 5, there is only a small region at large values of 
mi = (0.02 — 0.05) eV. Extending the analysis in Ref. |9], this is explained by the fact 
that the wash-out parameter Ki^. turns out to be 

K,, ^ ^ (^2 - mi) _ ^ ^^^^ ^ 

6 [2m2 + mi) 

while in the NO case Ki^. was proportional to mgoi. This constrains mi to be as 

large as possible, thus ruling out small values of mi. 

It is interesting to notice that in this case the allowed values for 6*23 lie in the sec- 
ond octant and correspond to the largest ones compatible with the current experimen- 
tal limits. The allowed values of the effective neutrino mass fall in a narrow range, 
mee = (0.05 — 0.07) eV. Therefore, 10 will be in any case fully tested from cosmology 
and /3/30z/ experiments during next years. We will see that this conclusion will hold also 
allowing Vl 7^ /• As usual, in the plots the red star corresponds to the minimum value 
of ^2 for which we have found a solution, 0:2 = 4.65. The corresponding set of values 
indicates approximately where the asymmetry has a maximum for a fixed 0:2 value. 
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Figure 4: Case Vl = I, 10. Scatter plot of points in the parameter space that satisfy the 
condition tib > 5.9 x 10^^ for a2 = 5 (yellow circles) and a2 = 4.65 (red star). In the top 
left panel the lower bound on Tjih (cf. eq. (|30|) ) is also indicated for the same values of 
^2 but with different symbols: 0:2 = 5 (grey squares), 0^2 = 4.65 (blue star). The three 
mixing angles are in degrees, the three phases in radiant. 
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Figure 5: Case Vl = I, 10. Plot of all relevant quantities versus mi for the set of values 
(^13 = 0.32°,^23 = 52.O3°,0i2 = 32.16°, p = 3.16,a = 3.48,5 = 2.47) corresponding to the 
red star in the previous figure (0:2 = 4.65). 

For this choice of values, in figure [5], we show the plots of the RH neutrino masses, of 
the asjTiimetry rjB, of the CP asymmetries 62, 620, of Ki, Kia and of K2, K2a versus mi. 
One can see how the heaviest RH neutrino mass M3 decreases with mi much faster and 
at mi ^ 0.001 eV one has M3 ~ 10^'' GeV. Therefore, as one can see from the central top 
panel, the CP asymmetries are this time strongly suppressed at mi 10^^ eV. On the 
other hand, in the range mi ^ (0.02 — 0.05) eV the CP asymmetries are large enough that 
successful leptogenesis is still possible. Notice, that this kind of solutions are a sort of 
modification of the solution obtained at large mi values for NO, simply shifted at somehow 
larger values. The asymmetry is therefore strongly depending on the initial conditions 
{K2T — 1). The first kind of solution, at small mi values, is completely absent. 

Therefore, though 10 is strongly disfavoured, it is not completely ruled out, a conclu- 
sion somehow very similar to that one obtained from completely independent arguments 
[26] . In this case, however, leptogenesis provides quite a precise quantitative test. 

We can conclude this section saying that these results confirm and complete those 
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shown in [9]. In particular it is confirmed that there are viable solutions corresponding 
to the different points shown in the figures falling in the currently experimentally allowed 
ranges of the parameters,. The model is therefore not ruled out. A further step is now 
to understand whether the model is predictive, excluding regions of the parameter space 
that future experiments can test. From the figures, as we have discussed, it is clear 
that assuming Vl = I such excluded regions exist and therefore one obtains interesting 
constraints. However, it is important to go beyond the simple condition = / in order 
to test the stability of the constraints for variations of V^. This is the main objective of 
the next sections. 



4 The case Vl = Vckm 

We now study how the constraints change when a misalignment between the physical 
basis where m/) is diagonal and the flavour basis, where the charged lepton mass matrix is 
diagonal, is considered, corresponding to Vl 7^ /. Since Vl is unitary, we can parameterize 
it similarly to the leptonic mixing matrix introducing three mixing angles, one Dirac-like 
phase and two Majorana-like phases. 



Vr 



' '-12 "-IS ^12*^13 ^13^ ^ 

qL r^L r^L ^16]^ „Z/ „Z/ ^16]^ qL 

■^12 "-23 '-12 *'23 *'13 '-12 '-23 "^12 ''23 ''13 ''23 '-13 

Vr,L r^L ^^L r^L ^i(5^ , 
■'12 ■'23 '-12 '-23 * 13 '-12 ■'23 ■'12 '-23 ■'13 '-23 '-13 / 



■diag {e'P^l^e'"'^) 



(41) 

where we defined sf. = sinOj"- and cfl- = cos 6'/"-. Therefore, we have now six additional 
parameters that give much more freedom. We will not explore the full parameter space 
but, in the spirit of S'0( 10)-inspired models, we will allow only small mixing angles at 
the level of the mixing angles in the CKM matrix. 

As a first definite example we repeat the analysis performed for the case Vl = I for 
a definite case where the are exactly equal to the mixing angles in the CKM matrix 
and therefore we set 9^^ = 0.21°, 6*23 = 2.3°, = 13°, where the latter is the measured 
value of the Cabibbo angle. 



4.1 Normal ordering 

For NO the results are shown in Figure 6. There is a first result to highlight: 0:2 values 
as low as ^2 = 1 are now allowed. This is an interesting result in connection with the 
study of realistic 5*0(10) models. At the same time this result also implies slightly lower 
values of M2 and consequentially of the minimum value of Tj-eh that can be now as low as 
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Figure 6: Case Vl = Vckm, NO. Scatter plot of points in the parameter space that satisfy 
the condition 77^ > 5.9 x 10"^ for a2 = 5 (yellow circles), a2 = 4 (green squares) and 
a2 = I (red stars). In the top- left panel, the same convention of Fig. 2 is adopted to 
indicate Trh- 
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^ 10^° GeV, as it can be noticed in the top-left panel in Fig. 6. In this case we have more 
generally calculated the minimum reheat temperature as 

r Mo Mo " 

T^/?~min , \ , , \ , (42) 

considering that in the case the asymmetry at the production can be either tauon dom- 
inated or e + /i dominated. This is because the third kind of solution that was highly 
suppressed in the case Vl = I, the right panel in Fig. 1, becomes now viable and is e + /x 
dominated, as we will discuss soon in more detail. 

Notice that if we compare the allowed points for eta = 4 with those found for Vl = I, 
the constraints on the low energy neutrino parameters are now less stringent. In particular 
an allowed region for values mi ~ 0.003 eV is also found for very small values of 6*13. 
Indeed, in the case Vl = I, and for small values of mi, the suppression of the wash 
out value Kir imposed a lower bound on 6*13. By choosing Vl = Vckm introduces the 
possibility of getting vanishing Ki^- even for zero 6'i3 angles. Extending the analysis of Ref. 
[9], one finds indeed that one configuration where Kij. is smaller than unity is attained if 
p = (mod 27r) and cos a = -[1/(12^^2)] Koi/K^a)] ~ -(5/12)(msoi/matm) ~ -10"^ 
This implies cr ~ 7r/2 (mod 27r), as confirmed by our numerical results. Including in the 
analysis the atmospheric neutrino mixing angle, as one can see from the panel with the 
constraints in the 6'i3 — ^23 plane, only values O2Z ^ 48° for 61^ < 10° are allowed for aa < 4. 
Notice that, for aa < 4, the allowed region in mi, 6*13, 6*23 only marginally overlaps, at small 
values of 6'23, with the region for the case Vl = I- This means that a measurement of 
these three quantities can distinguish between the two cases, Vl = I and Vl = Vckm, 
and not all combinations of these three quantities seem to be possible. We will be back 
on this point in the next section. 

In Figure 7 we plotted the relevant quantities for three particular choices of the pa- 
rameters, as indicated in the figure caption, corresponding to the three kinds of solutions 
found for Vl = Vckm- These three sets of values correspond to the three kinds of solu- 
tions that are found in the scan plots. The first two sets, corresponding to the left and 
central panels, give a tauon dominated asymmetry, while the third set, corresponding to 
the right panels, yields a muon dominated asymmetry. Notice that these three kinds of 
solutions are the same three kinds, with slight modifications, found for the case Vl = I- 
However, one can see that this time the third kind of solution, where the final asymmetry 
is muon dominated, also yields successful leptogenesis. The major difference that explains 
this result, is that for Vl = Vckm the flavoured CP asymmetries are not as hierar- 
chical as in the case Vl = I, as it can be clearly seen in the three panels showing the CP 
asymmetries in Figure 7. 
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Figure 7: Case = Vckm-, NO. Plots of the relevant quantities for three choices of the 
involved parameters. The long-dashed red lines correspond to a = r, the dashed blue 
lines to a = /i and the short-dashed dark yellow lines to a = e. Left panels: a2 = 4, 
= 1.7°, = 33.6°, 023 = 41.8°, 5 = 2.84, p =1.53 a = 3.24, pl = 0.12, = 2.56; 
central panels: aa = 5, ^13 = 3.3°, 6^ = 35.6°, 623 = 40.4°, 6 = -1.06, p = 2.87, a = 6.0, 
PL = 3.13, (XL = 3.25; right panels: 02 = 4, = 4.7°, 612 = 35.9°, 623 = 40.3°, 5 = -1.89, 
p = 0.065, a = 4.85, = 5.89, = 3.69. 




Figure 8: Case Vl = Vckm, NO. Constraints on the mixing angles without making use 
of the current experimental information from neutrino oscillation experiments. 

We have also repeated, as for the case Vl = I, the exercise to leave the mixing angles 
completely free, without imposing any experimental constraint finding the results shown 
in Fig. [HI One can see that in this case the points found when the current experimental 
constraints are imposed (the green points) fall in more marginally allowed regions, also for 
6*13. This might suggest that Vl = I seems to be a more attractive case than Vl = Vckm- 

4.2 Inverted ordering 

Finally, we also present in Figure 9 the constraints obtained for 10. Even though there is 
again a remarkable suppression of the allowed regions compared to NO, they are somehow 
less restrictive than for Vl = I- In particular now a broader range of values for mi is 
allowed and ^23 can be as low as ~ 45° for 0:2 < 4. This is also confirmed by the fact 
that lowest allowed value is now 02 = 2, much lower than in the case Vl = I (it was 
^2 = 4.65). However, it is still fair to say that the 10 case is only marginally allowed and 
certainly disfavoured compared to the NO case. 

5 Global scans 

The two specific cases that we discussed, Vl = I and Vl = Vckm, suggest an interesting 
sensitivity of S'0( 10) -inspired leptogenesis to slight deviations of Vl from the identity. 
This sensitivity was absent in the results found in A^i-dominated leptogenesis [5]. In this 
way it seems that one could even gain some information on Vl from low energy neutrino 
experiments. However, there is a potentially dangerous aspect of such a sensitivity: if 
for a slight variation of Vl the entire space of low energy neutrino parameters becomes 
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Figure 9: Case Vl = Vckm, 10. Scatter plot of points in the parameter space that satisfy 
the condition rjB > 5.9 x 10~^ for a2 = 5 (yellow circles), 0^2 = 4 (green squares) and 
02 = 2 (red star). 
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accessible, then any chance to test 5'O(10)-inspired leptogenesis is lost. On the other 
hand, from a comparison of the results obtained for the two definite cases, Vl — I and 
Vl — VcKM, one can understand that this does not happen. 

One can still suspect that for a continuous variation of the parameters in Vl, such 
that Vl changes from Vl — I to Vl = Vckm, new solutions appear so that any point in 
the space of the low energy neutrino parameters can be obtained for a proper choice of 
Vl. 

In this section we study this issue. We perform a global continuous scan of the pa- 
rameters for Vl between Vl = I and Vl = Vckm- Obviously a precise limit Vl = Vckm 
for such a global scan is somehow arbitrary. It should be therefore taken as a working as- 
sumption defining S'O(10)-inspircd leptogenesis, even more than the condition = 0{1) 
that, as we stressed many times, should not be regarded as a very restrictive assumption. 
Clearly within well defined realistic 5*0(10) models, more specific conditions on Vl should 
be obtained. In any case one expects that if the Vl satisfies the condition / < Vl < Vckm, 
then the allowed values for the low energy parameters should fall in the allowed regions 
for 5'O(10)-inspired leptogenesis. 

Therefore, in this Section we present the constraints on the low energy neutrino pa- 
rameters for a continuous variation of the values of the mixing angles in the range 
< < 9f^^ (i.e. for I <Vl < Vckm)- More exphcitly the shown scatter plots are 
obtained for the low energy neutrino parameters scanned over exactly the same ranges as 
for the case Vl — I. The three angles in Vl are scanned over the ranges < < 0.2°, 
< ^fa < 2.5°, < ^ 13°, while the three phases are scanned over [0, 27r]. In order to 
determine the allowed regions, we have followed the same strategy as in the case Vl — I, 
with a similar total number of scanned points, O(IO^). 

5.1 Normal ordering 

The results for NO are shown in figure 10. One can see how the allowed regions are 
approximately given by a super-position of those found for Vl = I and Vl = Vckm phis 
all intermediate solutions. The result is that now the correlations among the parameters 
found in the two special cases seem to disappear. There arc however still interesting non 
trivial constraints. What clearly survives is that the allowed points still cluster within two 
distinguished ranges of values for mi, one range at small values, nii ^ (1 — 5) x 10^^ eV, 
and one range at high values, mi ^ 0.03 — 0.1 eV, a distinction that is sharp for 0:2 = 4 
(green squares) while it is softer for 0:2 = 5 (yellow circles) . At the same time one can see 
that a global scan actually shows a slight correlation between mi and ^13 in the low mi 
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Figure 10: Global scan, NO. Scatter plot of points in the parameter space that satisfy 
successful leptogenesis {rjB > 5.9 x 10~^), for 0:2 = 5 (yellow circles), 0:2 = 4 (green 
squares) and 0:2 = 1 (red stars). 
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Figure 11: Global scan, NO, mi < 0.01 eV. Scatter plot of points in the parameter space 
that satisfy successful leptogenesis {rjB > 5.9 x 10~^), for 0:2 = 5 (yellow circles), 02 = 4 
(green squares) and 02 = 3.7 (red stars). The region below the dashed line in the left 
panel corresponds to the condition eq. P3|) . 

range while the interesting linear dependence between ^13 and 623 found for Vl = I seems 
now to be lost. 

However, it should be considered that these plots are projections on two-parameters 
planes of an allowed region in a seven-parameter space. Therefore, only a full multi- 
parameters analysis would be able to unreveal correlations involving more than two pa- 
rameters. Nevertheless, thanks to the distinct analysis that we carried out for the two 
special cases Vl = I and Vl = Vckm, one can catch sight of an interesting correlation 
among mi, 612, 613 and meg. To this extent, this time we have also plotted the constraints 
in the plane 613 — rriee, showing how the lower bound on rriee increases with 6*13. 

5.1.1 Low nil range 

In order to find out whether the linear dependence between 6*13 and 6*23 found for Vl = I 
(cf. eq. flMl) ) still holds for a global scan, we show in Fig. 11 the same constraints as in 
Fig. 10 imposing the condition mi < 0.01 eV, since the linear dependence was found in 
that range of values. We only show the constraints on the relevant parameters, therefore 
only those in the plane 613 — 623, in the plane 612 — O13 and in the plane ^13 — m^e- This 
time we could also easily find points for 02 = 3.7 (red stars), showing again how allowing 
for a Vl 7^ J the allowed regions get larger. 

One can see that the quite clear linear dependence eq. (l34l) between 6^13 and 6^23 holding 
for Vl = 1 1 now turns more, for the red star points at 0:2 = 3.7, into an allowed region 
below the dashed line showed in the figure and corresponding approximately to 



923 < 49° + 0.65 {9i3 - 5°) 
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(43) 



This result should be also understood in terms of the condition Kir ^ 1 (cf. (133|) ) when 
a very small Vl is allowed clearly yielding a dispersion around the linear dependence 
eq. (Elj). Notice that inside this region there are still sort of sub- regions that seem to be 
excluded. 

We can summarize these results saying that, at low values of rrii < 0.01 eV, there is 
an interesting testable constraints in the plane 613 — 623 given by the relation eq. fH3|) . In 
particular experiments that are already taking data such as the nuclear reactor experiment 
DOUBLE CHOOZ [27j and the long baseline experiment T2K [2B] have the capability of 
a 3cr discovery of values 6*13 > 8°. Our results seem to suggest that if such high ^13 values 
will not be found, then a restricted range of values for 6*23 is predicted. For example, if 
^13 < 8° then 623 < 51°, and if ^13 ~ 6° then 623 < 48°. Such a constraint on 623 should 
be also tested during next years with quite a good accuracy by the T2K experiment [28] . 
These constraints in the plane 6^13 — 6*23 should be considered at this level indicative, and 
should also consider that they are quite sensitive to the value of 0^2. 

Notice that at the same time, cosmological observations and/or neutrinoless double 
beta decay experiments should also be able to test the condition mi < 0.01 eV. It should 
be therefore appreciated that this scenario will be tested during next years. 

It is also interesting to notice (see right panel in Fig. 11) that there is a linear depen- 
dence between and 6*13 as well. In particular, for 02 < 4, at large values 613 > 6° 
one has rriee — 10~^eV and even for 613 > 8° one has rriee — 3 x 10^'^ eV. These values 
for rriee are below the sensitivity of future planned experiments (> 0.01 eV) such as EXO 
[29] . However, at least, rriee cannot be arbitrary small but has a lower bound that, for 
sufficiently large 613 values, is 3 times below the currently planned reachable experimental 
sensitivity, a very small value but maybe not completely hopeless. 

Within the two-parameter analysis we are presenting, we cannot draw sharper pre- 
dictions but is seems quite plausible that from a more involved multi-parameter analysis 
precise correlations could emerge, maybe also involving the solar neutrino angle 612- In 
this respect the central panel in figure 11 suggests that the solar mixing angle could indeed 
play also a role and that maybe sharper predictions in the 3 parameter space (^^13, Ou, 623) 
exist. 

5.1.2 Large mi range 

We can also study how the allowed regions would reduce requiring large values mi > 
0.01 eV. The results are shown in figure 12. One can see that in this case one obtains 
very clear constraints that will allow to test this scenario during next years in a quite 
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Figure 12: Global scan, NO, mi > 0.01 eV. Scatter plot of points in the parameter space 
that satisfy successful leptogenesis {rjB > 5.9 x 10~^), for 0:2 = 5 (yellow circles), 0:2 = 4 
(green squares) and 0^2 = 1 (red stars). 

unambiguous way. First of all from the Fig. 10, thanks to the very precise values of the 
Majorana phases, one can notice that there is a very clear relation between mi and rriee- 
Second, one can see from the left panel of Fig. 12 how there is an upper bound 623 < 46° 
for a2 < 4. For values of ^13 ~ (5 — 6)°, one has even 623 < 41°. It should be said however 
that at these large mi values, one typically obtains a final asymmetry that depends on the 
initial conditions. Since we are assuming vanishing initial N2 abundance and vanishing 
initial asymmetry, these constraints should be regarded as the most stringent ones, but 
likely also the best motivated ones. 

5.2 Inverted ordering 

Finally, we repeated the global scan for 10 as well and the results are shown in Fig. 13. 
One can see how the allowed regions somehow merge those found for the two extreme 
cases Vl = I and Vl = Vckm- There is therefore nothing really new. 10 is quite strongly 
constrained and it will be fully tested in next years. In particular we can notice again 
how there is a clear lower bound on 623 rather than an upper bound as in NO. More 
particularly, one can notice that the allowed region in the plane mi — 623 is approximately 
described by 

023 ~ 43° + 12° log(0.2 eV/mi) (44) 

(the dashed line in the upper right panel). It is then quite interesting that S'O(IO) inspired 
leptogenesis is able to distinguish NO and 10 even at mi > 0.01 eV, when the same values 
of mee and of ^ • mj (the quantity tested by cosmological observations) are found both 
for 10 and for NO. From this point of view 5*0(10) inspired leptogenesis provides a way 
to solve this ambiguity. 
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Figure 13: Global scan, 10. Scatter plot of points in the parameter space that satisfy 
successful leptogenesis {rjB > 5.9 x 10~^), for 0:2 = 5 (yellow circles), 0^2 = 4 (green 
squares) and 0:2 = 1.5 (red stars). The dashed line is the eq. (jSj). 
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6 Final remarks 



We have derived constraints on the low energy neutrino parameters from S'O(10)-inspired 
leptogenesis. Our investigation shows that even minimal leptogenesis, based on a type I 
seesaw mechanism and assuming a thermal production of the RH neutrinos and with a 
traditional high mass scale RH neutrino spectrum, can be testable within a well motivated 
framework, where the see-saw parameter space is restricted by the 5'O(10)-inspired con- 
ditions. The role played by the N2 decays is crucial in this respect, not only in re-opening 
the viability of these models. The presence in the A''2-dominated regime of a double stage, 
a production stage and a lightest RH neutrino wash-out stage, seems to introduce, as 
shown simultaneously both by the numerical and by the analytical results, a strong di- 
rect dependence on neutrino mixing angles as well, in addition to the dependence on the 
absolute neutrino mass scale, already found in usual A''i-dominated leptogenesis [2]. 

Interesting predictions, that can be tested in future years, with intriguing correlations 
involving the absolute neutrino mass scale and the neutrino mixing angles emerge. 

In the significant case of NO with low uii values, the neutrinoless double beta decay 
effective mass seems to be too small to be measured but not arbitrary small and in any 
case future experimental results can be anyway useful to restrict the allowed regions for 
the other parameters and sharpening the predictions. 

The results for Vl 7^ / seems also to be sensitive to Vl itself and they therefore suggest 
that there is an opportunity to gain information on it, an interesting point within studies 
of specific 5*0(10) models. It is quite interesting that there is an allowed region in the 
parameter space that allows large values of 6^13 testable with on-going reactor neutrino 
experiments and that for these large values the models favours either large or small 623 
values depending whether mi < 0.01 eV or mi > 0.01 eV. 

In the small mi range it is also interesting that the constraints are completely inde- 
pendent of any assumption on the initial conditions, a point that maybe makes this option 
more attractive. It is actually quite interesting that this conclusion is also supported by 
completely independent and general considerations based on the possibility to reproduce, 
without a particularly fine tuned Ur matrix, the observed atmospheric to solar neutrino 
mass ratio, matm/''^soi — 6, starting from hierarchical neutrino Yukawa couplings. It is 
found [30] that this experimental observation is far more natural if the lightest neutrino 
presents a much stronger hierarchy than the the two heavy ones, as it occurs in the region 
that we have found at small mi. It should be also stressed again, that since our results 
are independent of 0:3 and ai, as far as M3 > 10^^ GeV and Mi < lO^GeV, they hold 
even for a Yukawa couplings hierarchy milder than in the case of up quark masses. This 
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can help to make even more natural to reproduce the result matm/''^soi — 6 without a fine 
tuned Ur. 

A more precise measurement of 612 could also play a relevant role in testing these 
models, a point that should be addressed by a more involved multi-parameter analysis. 
A future accurate determination of the neutrino mixing angles will be therefore crucial 
to test S'O(10)-inspired leptogenesis and could even yield some interesting information 
on the matrix Vl. In conclusion, it seems that S'O(10)-inspired leptogenesis provides an 
interesting well justified example that gives some hopes about the possibility of testing 
minimal leptogenesis even only with low energy neutrino experiments. It will be then 
quite interesting in next years to compare the experimental results with the constraints 
and the predictions from 5'O(10)-inspired models that we discussed. 
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